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Abstract. We exam various geometric active contour methods for radar im- 
age segmentation. Due to special properties of radar images, we propose our 
new model based on modified Chan-Vese functional. Our method is efficient 
in separating non-meteorological noises from meteorological images. 



1. Radar image processing 

Weather radar data quality control is extremely important for meteorological 
and hydrological applications. For weather radars, scatterers in the atmosphere 
are not only meteorological particles like cloud, rain drops, snowflakes, and hails, 
but also non-meteorological particles such as chaff, insects, and birds. For radar 
meteorologists it is a major issue and challenge to design numerical scheme which 
can automatically and subjectively distinguish meteorological echoes from non- 
meteorological echoes. Non-meteorological echoes can contaminate radar reflectiv- 
ity and Doppler velocity measurements, and subsequently cause errors and uncer- 
tainty on radar data applications in quantitative precipitation estimation, as well as 
in assimilation in numerical model for weather prediction. Automatic detection of 
tornado or mesocyclone vortex among meteorological echoes is certainly another big 
challenge and has great potential in improving severe weather forecast and saving 
human life. 




(a) The famous May 03, 1999 Tornado (b) Noised storm image: radar noise is em- 
(Moore, Oklahoma) radar image bedded in storm image 

Figure 1. Challenges in radar images 
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For different echoes, tlieir graphic properties such as pattern, intensity and tex- 
ture, are different. See, for example, fig. [H Such differences enable us to design 
new active contour model to automatically extract the most distinguishing graphic 
properties from echoes, and to segment them easily from other noises. Our methods 
shall also be very useful in automatic detection of tornadic supercell, as well as in 
storm classifications and tracking in the future study. 



2. Active contour 

Active contour is the procedure that we deform a given curve so that a given func- 
tional of the curve will achieve its local minimal value. This method is wildly used 
recently in computer vision in seeking the edges or contours of given images. See, 
for example, Mumford and Shah [4], Kass, Witkin and Terzopoulos [3], Caselles, 
Kimmel and Sapiro p], and Chan and Vese [2]. 

Let uq{x, y) : — > i? be the gray level function of a given image. If uq is smooth, 
then the edge of the image are those points {x,y) where |Vuo| is relatively large. 

The geometric contour and snake models aim to detect edge automatically based 
on the size of |Vito|. Let C(s) : [0, L] — > be a closed curve, where s is its arc 
length parameter. One can introduce an edge-detector function g : Q ^ so that 
g{z) — > as 2 — > oo. A typical example of such function is given by 

We define the energy functional of C by 

(2.1) h{C) f\{\7uo{C{s)))ds, 

Jo 

then to find the edge of image uq can be reduced to seek the local minimal for Ii 
(the geometric active contour model '1 ): 

(2.2) h{edge)^infch{C). 

The snake model 3 is to introduce, for a parameterized curve C(p) : [0, 1] ^ M^, 
the following energy functional : 

(2.3) l2{C)^af \C'{p)\^dp + /3 f \C"{p)\dp^\f \Vuo{C{p))\^dp, 

Jq Jq Jo 

where a, (3, X are all positive parameters. The first two terms represent the internal 
energy of the image, which usually are used to smooth the curve; The third term 
represents the external energy, serving as the indicator for edge. The edge of the 
image then can be found by minimizing l2'. 

(2.4) /2 (edge) =inf /2(C). 



To automatically detect the edge via an iteration scheme, one introduces a family 
of curves C{p,t) : [0, 1] x [0, oo) and the deformation path. For example, for 

active contour model (12. 2p the curve evolution (gradient flow equation) is given by 

(2.5) Ct = {kg - Vg ■ N)N, 

where k is the curvature function and N is the inner unit normal vector of curve 
Cip,t). 
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Numerically, such iteration can be realized via the powerful level set method of 
Osher and Sethian [5]. Embed C{p, t) as a nodal line of a smooth function ^{x, y, t): 
C = {(x, y, t) : $(a;, y, t) = 0}. From 9f $((7, t) ~ 0, we are led to evolve $ by 

— =5(|V.o|)d*-(^)|V$|+<V5,Vci>> 

<^{x,y,0) = $o(a;,y), 

where '^^{x^y) is the initial level set function. In practice, one can choose $0(2;,?/) 
to be a signed distance function to a given initial curve C(p, 0). Fig. 2 (b), fig. 3 
(b), and fig. 4 (b) show the result for image processing based on such scheme. 

If the given image uo(x, y) is not smooth, the edge of the image is not well defined 
based on the derivative of the gray level function. The human being's perspective 
for the edge of a non smooth image basically is to identify the boundary of different 
groups. To identify such boundary, one can use Chan-Vese energy ^2j: 

/3(C,ci,C2) : = / \uq - ci\'^ dxdy + / \uf) - 02]"^ dxdy 

J inside{C) J outside{C) 

(2.7) + /I ■ {length{C)) + v ■ {Area{inside{C))) , 

where ci, C2 are constants to be adjusted in iteration, /i and v are fixed parameters. 
The last two terms are smoothing terms. The edge is again sought by minimizing 
/3(C, ci,C2): 

(2.8) h{edge,ci,^,C2,*) ^ inf h{C , ci , C2) . 

C,Ci ,C2 

Again, numerically level set method can be used for such deformation. Introduce 
the Heaviside function and its derivative 

Hi^)^{ J' 'I'll Si.)^±Hi.). 
[ 0, if z > 0, dz 

Embedding C{p, t) as a nodal hne of a smooth function <I>(a;, y,t): C — {{x, y, t) : 
$(a;, y, t) — 0}, we can re- write the energy functional I^{C, ci, C2) as 



J3($,ci,C2) : = / \uo-ciYH{<^{x,y))dxdy+ / \uo - C2WI - H{^{x,y)))dxdy 
(2.9) +n S{<^{x,y))\W<^{x,y)\dxdy + iy / H{<^{x,y))dxdy. 



For fixed $, minimizing Js^^, ci, C2) with respect to c,; yields 

ci(<i>) = average(uo) in {$ < 0} 
C2(<i>) = average(uo) in {$ > 0}. 

Once ci and C2 are fixed, we minimize J3 via deforming $ along the gradient 
direction: 

-^=6i^){^i-div{^^)^u~{uo-c^m' + {uo-C2m^} in (O,oo)xl], 
$(a;,?/,0) = $o(a;,y) in n, 
— = on ail 
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where ^o{x, y) is a signed distance function to a given initial curve C(p, 0). Fig. 2 
(c), fig. 3 (c), and fig. 4 (c) show the result for image processing based on such 
scheme. 



3. Modified model 

As being pointed out in [2] , The Chan-Vese model is originated in the Mumform- 
Shah model [4j: 

C)^ ^1- length{C) + A / \uo{x, y) - u{x, y)\^dxdy 

\Vu{x,y)\^dxdy, 

n\c 

where mq : — > M is the given image, K C is the curve that we will deform. The 
sharp boundary of uq, where |Vuo| is large or discontinuous can be detected via 
minimizing F^'^^ {u^C). Roughly speaking, if we replace m by a constant function 
in the above, we can see the prototype model similar to that of Chan-Vese. These 
models are all viewed as minimal partition problems. 

We observe that radar noises usually have relatively low intensity to severe 
storms, and radar signals for storm are usually uniform in certain region. To seg- 
ment the severe storm from usual radar noises, we introduce the following modified 
Chan-Vese functional: 

/4(C, c) : = \uo — a ■ M\^dxdy + \uo — c\^dxdy 

J insideiC) J outside(C) 

(3.1) + /I ■ {length{C)) + v ■ {Area{inside{C))) , 

where M — maxajgo Mo(a^), c is a constant to be adjusted in iteration, a, /i and i/ 
are fixed parameters. In practice, parameter a can be determined by comparing 
the maximal intensity near radar and the maximal intensity of storm. 

Embedding C{p, t) as a nodal line of a smooth function $(a;, y, t): C ~ {{x, y, t) : 
^{x, y, t) — 0}, we can re- write the energy functional Ii{C^ c) as 

J4($,c):= / \uq- a- M\'^H{^{x,y))dxdy + i \uq - c\^{l - H{^{x,y)))dxdy 
Jn Jn 



(3.2) +fi / Smx,y))\V^x,y)\dxdy + iy / iJ($(a;, y))dxdj/. 

Jn Jn 

For fixed $, minimizing Ji{^, c) with respect to c yields 

c($) average(uo) in {$ > 0}. 

To derive the first variation of the functional, we consider slightly regularized 
version of functions i?e and H'^ — such that G C^fl), ^ H and 5, 
and the modified functional: 

J4,e($,c):= / \uo-a-MfH,mx,y))dxdy+ [ \uo - c{<^>)f{l - H,mx,y)))dxdy 
Jn Jn 

(3.3) 

+ /X / d,{'^{x,y))\S/'^{x,y)\dxdy + ly / H,(^{x,y))dxdy. 
Jn Jn 
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Its first variation is 

<(5J4,e,V>= / iuo-a-MfS^{(i>{x,y))'iJjdxdy- / {uo ~ c{^))^S,{<!>{x,y))'ipdxdy 
Jn Jn 

+ /i / — -Vijdxdy + i^ / S^mx,y))iljdxdy. 

Jn |V$(a;,y)| 

Thus we can minimize J4 c via deforming $ along its gradient direction: 
(3.4)^^ 'v$ 

— = S,{^){^^ ■ divi-^^) -iy-{uo^a- Mf + {u^ ~ c{<^)f] in (0, 00) x 

$(x,y,0) = $o(a;,y) in f7, 

on oil. 



iv$i a? 



4. Experimental Results and Discussion 

An image of reflectivity factor field of tornadic supercells, observed by a S-band 
weather radar near Oklahoma City, Oklahoma, is presented in Fig. 1(a). The violent 
tornadoes generated from the weather system ripped through Oklahoma and Kansas 
and killed 48 people while demolishing houses and business. It had caused at least 
$500 M in property damage (http:/ /www. nssl.noaa.gov/headlines/outbreak.shtml). 

We first use geodesic active contour model (|2.6p for the image. It results in 
catching all boundaries, including radar noises (fig. 2 (b)). We then apply standard 
Chan-Vese model (I2.10p with /i = 5, i/ = (fig. 2(c)). Still it keeps almost all 
boundaries from radar noises. We finally apply our model (|3.4p with a — 0.7 in fig. 
2 (d). The result shows that almost all radar noises are successfully skipped. 

Fig. 3 (a) usually is a challenge radar image for processing. The radar noises are 
embedded in storm image (in fact, radar is underneath the cloud). Geodesic active 
contour is very sensitive to the initial curve. It usually contracts curve. Fig. 3 (b) 
is a failure via geodesic active contour. Fig. 3 (c) is the result using Chan-Vese 
model (choose — b, v = Q); Fig. 3 (d) is based on our model p.4p (with a — 0.4). 
There is no big difference between fig. 3 (c) and fig. 3 (d). 

Fig. 4 (a) is another radar image with radar noise separated from storm. Chan- 
Vese can not remove radar noise completely. Our model with suitable parameter ( 
a — .6) works fine. 

Finally, we compare the results using Chan-Vese model and our model with 
different parameters. 

First we consider Chan-Vese model with different parameters: 

(4 1) I ^^5m{^i■d^v{^^)~v~X{u,-c,mf + {u,-c,mf] 

y $(x,?/,o) = 0o(x,?/), 

where A is a positive parameter. The results using Chan-Vese model with different 
A are presented in fig. 5. 

Next we compare the results using our model with different a in fig. 6. It can 
be seen that for a in certain range, our results are relatively stable. Therefore, for 
different a, if we let 

c°"*($) = average(uo) in {$ < 0}, 
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we can develop a program which can automaticaUy determine which a we shall 
choose based on the changing of c°"*(<i>) as a changes. 

5. Conclusion 

We compare various models and their applications to the segmentation of radar 
images. We propose our new model. Our method is more efficient in outlining more 
severe storm images, and skipping the usual radar noises. 

ACKNOWLEDGMENT. M. Zhu is partially supported by the NSF grant 
DMS-0604169. 



Iterations 2000 Iterations 




(a) Original Tornado image (b) Geodesic active contour: failed 



2000 Iterations 2000 Iterations 




(c) Chan-Vese model: Catch almost all (d) Our model: Catch serious storm, skip 
boundaries, majority of radar noise is kept radar noises 

Figure 2. Comparison for different method 
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(a) Radar noises are embedded in storm (b) Geodesic active contour: Can not ex- 
image pand 



2000 Iterations 20O0 Iterations 




(c) Chan-Vese model: skip radar noises (d) Our model: Catch serious storm, skip 

radar noise 



Figure 3. Comparison for different method 
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(a) Radar noises are away from the storm (b) Geodesic active contour: Catch all 

boundary 





(c) Chan-Vese model: Catch almost bound- (d) Our model: Catch serious storm, skip 
ary, small part of radar noise is kept radar noise 

Figure 4. Comparison for different method 





(a) Original Tornado image 



(b) A = 2 





(c) A = 3 

Figure 5. Chan-Vese 



(d) A = 4 



model with different A 




(c) a = 0.7 (d) a = 0.9 

Figure 6. Our model with different a 



